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ABSTRACT 

The addition of m singlet right-handed neutrinos to the Standard Model leads to radiatively 
generated mass corrections for the SU{2)l doublet neutrinos. For those neutrinos which 
are massless at the tree level after this addition, this implies a small mass generated at the 
two-loop level via W^ exchange. We calculate these mass corrections exactly by obtaining 
an analytic form for the general case of n doublets and m singlets. As a phenomenological 
application, we consider the m = 1 case and examine the masses and mixings of the doublet 
neutrinos which arise as a result of the two-loop correction in the light of experimental data 
from two sources which may shed light on the question of neutrino masses. These are (a) 
the neutrino detectors reporting a solar neutrino deficit (and its resolution via Mikheyev- 
Smirnov-Wolfenstein matter oscillations), and (b) the COBE satellite data on the non-zero 
angular variations of the cosmic microwave background temperature (and its possible impli- 
cations for hot dark matter). Within the framework of the extension considered here, which 
leaves the gauge group structure of the Standard Model intact, we show that it is possible for 
neutrinos to acquire small masses naturally, with values which are compatible with current 
theoretical bias and experimental data. 



1 Introduction 

It is fair to say that the problem of understanding the origin of fermion masses is one of 
the most perplexing questions facing particle physics today. The Standard model |I| can 
reproduce the observed fermion masses via electroweak symmetry breaking and the Higgs 
mechanism, but provides no explanation for their values. When such an understanding is 
obtained, one of the issues that it must clarify is the smallness of neutrino masses (if, indeed, 
neutrinos are massive) relative to those of the other fermions. An attractive explanation 
for this observed feature of the fermion mass spectrum is the see-saw mechanism 0. It 
postulates the existence of right-handed neutrinos with masses of the order of the next 
energy threshold and uses this in combination with the Higgs mechanism to generate light 
(Majorana) neutrino masses via an effective dimension five operator. 

Given our present ignorance of the origins of mass and the lack of experimental pointers 
towards any particular mechanism, it is important to keep an open mind on the smallness 
of neutrino mass. In this paper, we explore, via detailed calculation, the issue of radiatively 
generated neutrino masses, since this is also a natural way in which masses small compared 
to those of other fermions may be generated. 

Any such effort needs to invoke physics beyond the Standard model. In view of the 
extraordinary and demonstrated robustness of the model to experimental tests over the last 
twenty years, we have thought it reasonable to make the simplest possible extension to 
the standard theory and study its effect on neutrino masses via radiative corrections, i.e. 
the addition of m SU{2)l(^U{1)y singlet right-handed neutrinos. A priori, there is no 
connection between their number and that of the doublet neutrinos, hence the simplest case 
corresponds torn = 1, i.e. the addition of one right-handed singlet neutrino to the standard 
model §. 



The gauge group structure of the weak sector remains unchanged as a consequence of 
this extension, but majorana mass terms incorporating the scale of new physics are now 
allowed. We do not speculate on their origin, but only note that it would require invoking 
an additional global symmetry (such as a conserved lepton number) to set these to zero. 
The doublet neutrinos acquire radiative (and, in some cases, tree-level) masses due to the 
presence of the singlets, as we discuss below. The radiative masses arise (via mixing) due to 
a two-loop mechanism 0, ^ involving the exchange of W^ bosons. In Sections II and the 
Appendix, we calculate, exactly and in analytic form, the two-loop masses accquired by the 
initially massless doublet neutrinos. Our calculation is general and valid for any number of 
doublet and singlet neutrinos, but in order to obtain phenomenologically useful information, 
we focus, in Section III, on the m = 1 case. Even this simplest extension of the Standard 
model introduces four new parameters into the theory. On the issue of neutrino masses, it 
is non-accelerator experiments that provide information on the cutting edge. Hence we have 
chosen to examine the results for the 171=1 case in the light of (a) the MSW solution to 
the solar neutrino deficit seen by the Kamiokande [0 , GALLEX [^ , SAGE and Homestake 



[[10| neutrino detectors and (b) the implications for hot dark matter (neutrinos) from the 



recent COBE observations on the anisotropy of the microwave background [|TT|. Invoking 
this experimental information restricts the parameter space and consequently, in the context 
considered in this paper, permits a handle on the range of the mass scale characterizing 
physics beyond the standard model. We show that doublet neutrino masses compatible with 
both (a) and (b) above can result naturally from such physics at the several hundred GeV 
scale. 



2 Radiative Generation of Neutrino Masses 

In this section we give a description of an exact general procedure for calculating two-loop 
neutrino masses applicable to any extension of the Standard model which incorporates singlet 
right-handed neutrinos. (We remark below on the reason why a one-loop mass does not arise 
in the sitiuation considered here, where only right-handed handed neutrinos are added to the 
existing particle spectrum. ) After setting up the generic integral that needs to be calculated 
we describe the procedure for evaluating it exactly in the Appendix. 

The lepton sector of the extension considered here has, in general, n (> 3) doublet fields 
I'^'iL hbi'" and m singlet fields {v'alY = {'^a)r- ( Here i = l....n, A = l....'m and u'^ = Cu'^ 
is the charge conjugate spinor.) In addition, one has the charged lepton SU{2)l singlet 
fields liR. The primes on the neutrino fields denote weak eigenstates as opposed to physical 
particle states. Without any loss of generality, we have assumed that the weak eigenstates 
li are the same as the corresponding mass eigenstates, i.e. the charged lepton mass matrix 
is diagonal. 

As noted in the Introduction, in addition to the Dirac mass terms, the most general 
Lagrangian consistent with the gauge symmetry of the Standard model also contains possible 



Majorana mass terms for neutrinos of the form TnABi^'AhY^'BL- ^"^ ^^e minimal model under 
consideration here, such terms must be bare mass terms, but in a more involved model they 
could arise, for instance, due to the vacuum expectation value of a singlet higgs. To facilitate 
discussion, we combine all the left handed neutrinos into a (n + m)-dimensional vector in 
the flavour space denoting it by z/^^. The most general mass term is thus given by 

n n+m 



C^ = Y.l^i^^LliR+ E «LTMapu'f^L + h.C (2.1^ 

j=l Q,/3=l 



Here A^ is a complex symmetric El (ra + m) x (n + m) matrix of the form 

M = 



^nxn ^nxm 



f2.2l 



-^mxn ^'^mxm, 



with D and M denoting the Dirac and the Majorana mass terms respectively. The first 
block is identically zero in the absence of a non-trivial vacuum expectation value for a 
S'f/(2)/^-triplet higgs field. (This restriction is imposed not only by our philosophy of minimal 
extension, but more importantly, by my/ /mz — the observed ratio of the gauge boson 
masses.) M. can be diagonalized by a biunitary transformation of the form 



V'^MV = M = diagonal (mo.) 



The mass eigenstates ( Va) are then easily identified to be 



l^L 



vh 



(2.3) 



(2.4) 



The relevant piece of the weak Lagrangian is then given by 



>C„k = J^.w^ 



(2.5) 



where 



n n+m 



^; 






10 . . — 7 . . . ' — 



ig 



'^Cw i=i 



fL 



2cw 



n+m 



Y^ klfiPLk + Y {K^K)apUo,'yf,PL^p 

a,t3=l 



i=l 



(2.6) 



Here cw = cos6'vi/, where 9w is the Weinberg angle, g = e/sin6'v^, P^ = (1 — 75)/2 and 



K 



-^nxn U 







V 



(2.7) 



^That A4 has to be syminetric is evident from the charge conjugation property of fermion bihnears. 



is the (n + ?72)-dimensional analog of the quark sector Cabibbo-Kobayashi-Maskawa matrix. 
Note that though KK"^ = diag(/„xn5 0), {K'^K)^^ ^ 6a/3- Thus we do indeed have flavour 
changing neutral currents (FCNC) in the neutrino sector. 

Having set up the general formalism, let us now concentrate on the case where n > m. 
There exist then n — m neutrinos that are strictly massless at the tree level. We now calculate 
the changes to such a spectrum accruing from quantum corrections. 

Before proceeding, in view of the fact there exist massive neutrinos and also FCNC's in 
the neutrino sector, it is appropriate at this point to remark on the possibility of one-loop 
graphs with a Z or a Higgs exchange introducing a non-trivial correction to the neutrino mass 
matrix. However, it can be easily seen that it is possible to rotate the neutrino states such 
that only m of them have Yukawa couplings to the Higgs. Thus, only those doublet states 
that are massive at tree level obtain a Higgs induced mass at the one-loop level. In addition, 
since the flavor- changing Z couplings have the same mixing parameters as the flavor-changing 
Yukawa couplings, the one-loop Z exchange diagrams do not contribute to the masses of the 
n — m neutrinos which are massless at the tree level. This reasoning applies at all orders to 
any diagram where all virtual particles are neutral. Hence the relevant diagram to compute 
is that given in Fig.(l). 

We shall work in the weak interaction basis for the external neutrinos and the mass basis 
for all the virtual particles. Furthermore, we shall concentrate only on the first n x n block 
of Ai, i.e. on the generation of Majorana mass terms for the doublet neutrinos. In the 
unitary gauge, the correction to the neutrino propagator is then given by 

-<-^--P^ ^^J ,,7 ^-P^ (2-8) 






q^ — m^r k"^ — m^ 



w 



The mass correction is of course given by Mi/ = ^ij (p = 0), and after some algebra this 
leads to 



4 n+m 



^if = TE^"<-^i. 



a=l 



(Tk drq {k + qyk ■ q 



q'^ + k^ 



4 + 






-4 



m 



w 



(2.9) 



where 



%;« = ik + qf {{k + qf - ml) (q' - f,^){q' - ml){e - m^)(P - li]) 



(2.10) 



We see that the mass corrections would be identically zero if iria = 0, Va. This ought to be 
so as any mass renormalization must be proportional to the bare mass terms. The integral 
above has a naive degree of divergence of 4. However, note that 

n+m 



Y: m^K^Kl = M., = 



a=l 



and hence 



n+m 



Y: m^KlKl 



{k + qf 



n+m TV-t TV-t ™,3 
a 



(2.11) 



(2.12) 



a=i (k + qf - ml ^^ {k + g)2 - ml 

This clearly is analogous to the GIM mechanism in the quark sector. Even on substitution 

of eqn(^]T^) in eqn( p.9D , the integral in the latter is still formally divergent. Notice, however. 



that this is but an artifact of the unitary gauge and is not a real divergence |I2|. In fact, by 
invoking identities similar to eqn( |2.12| ) or equivalently, by working in the Feynman gauge, 
one obtains □ 



n+m 






t Z^t 



a=l 



4 + 4^i±^+# 



ra 



w 



m 



w 



A / 2 2 2 n 2 



w) 



(2.13) 



where 



A{m1, 1712, "^3? "^4? "^5? ^ 



6) 



^This has often been cited in the hterature as a GIM-hke cancellation, but in our view the two are quite 
different. 
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d'^k d'^q k ■ q 



{q^ + ini){q^ + rn^ {{k + qY + m^} {{k + qY + ml} {k^ + ■ml){k'^ + ml) 

(2.14) 

is an Euclidean integral evaluated in the Appendix. 

The expression in eqn( |2.13| ) thus represents the Majorana mass generated for the doublet 



neutrino at the two-loop level. In operator language, it arises from terms of the form 



(L,i)'=L,x005 (2.15) 

where Ln represent the doublet lepton fields, is the usual higgs field and S represents the 
lepton number violating operator (whether a singlet higgs or a bare mass term). We note 
that this five dimensional effective operator for the radiative masses is the same as that for 
the conventional see-saw mechanism. The difference between the two resides in the scale 
of mass generation. Two-loop radiative masses compatible with the solar and COBE data 
can arise from right-handed neutrinos at the several hundred GeV scale, as we show below, 
whereas the see-saw mechanism generates similar valued masses via heavy neutrinos at the 
grand unified scale. 

We also note that though the corrections ostensibly are proportional to m|^ (eqn.( |2.1^ )), 
the actual dependence is linear (apart from logarithmic corrections) due to suppressions 
hidden in A. As rria becomes larger and terms of the order of {fii/maY become negligible, 
the correction goes as T^aKaiK'^jiria, which is simply the {ijY^ element of the tree-level mass 
matrix, and hence zero for the cases of interest here. 



Finally, we remark that a complex A^ in eqn. ( |2.1| ) obviously leads to a complex digonal- 
izing matrix V and hence possibly to CP-violating processes. However, since there is no 
evidence as yet of any such non-conservation in the leptonic sector, we have, in the interests 
of simplicity, chosen to perform all numerical calculations assuming a real neutrino mass 
matrix. 
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3 Application: The Solar Neutrino Deficit and COBE 
Data 

In order to make a connection to experiment and phenomenology, we now specialize to the 
n = 3 and m = 1 case and examine the two loop mass corrections in the context of (a) the 
MSW solution to the solar neutrino deficit reported by various detectors ^ |, |^, |T0[ and 



(b) recent COBE [jTT| data and its implications for neutrinos as dark matter. 

The solar deficit is the only long-standing possible evidence for physics beyond the Stan- 
dard model, and the MSW mechanism is its most popular resolution. In its essence, the 
mechanism requires neutrinos to be massive (and non-degenerate), allowing the interaction 
eigenstate z/g (assumed to comprise predominantly of the lightest mass eigenstate) to oscillate 
to z/^ or z/t- due to the difference in the forward scattering potential seen by the two states 
in their passage through solar matter. It thus identifies a range of vacuum mixing angle and 
mass squared difference values which are compatible with the deficit observed by the various 
detectors. Figure 2, excluding curves labelled (a), (b) and (c), is taken from Ref. [1^ and 



shows the familiar two-fiavor mixing MSW solution space, where 6 is the Cabibbo mixing 
angle and Am^ is the difference of the squares of the two neutrino masses, which, in the 
present context, are acquired at the two- loop level. 

COBE data on the anisotropy of the microwave background, while not making a definitive 
statement on the nature of dark matter, seem to suggest that it may have both hot and cold 
components, with the former being a neutrino (since it is the only known hot dark matter 
candidate) with a mass of ~ 10 eV. 

We use both of the above considerations to restrict the rather large parameter space 
available to us. 

In the scenario with one additional singlet, we have two massive and two massless neu- 



trinos at tree-level. The two massive ones acquire both one-loop and two-loop corrections, 
which we neglect, and the massless states acquire small masses at the two-loop level. The 
two tree-level masses and all the radiative corrections are expressible in terms of four in- 
put mass parameters for the matrix Ai. For various plausible (fixed) values of rria, (the 
singlet mass, signifying the scale of new physics) and the added constraint that the other 
neutrino with a tree level mass lie in the 10 eV range, we obtain a one parameter set of 
curves (see Fig. 2) which denotes the intersection of the "two-loop space" with the MSW 
solution space. Note that restricting ourslves to the two dimensional MSW space imposes 
a third constraint, i.e. that the z/g mixes predominantly with only one other state. Curve 

(a) in Figure 2 corresponds to a singlet mass of 100 GeV and a u^ mass of ~ 8.6 eV. The 
two-loop masses and mixings of z/g and z/^ are then such that they span the MSW space as 
shown. Curve (b) corresponds to a singlet mass of 400 GeV and a z/^ mass oi ^ 7 eV. z/^ 
and z/e then acquire radiative masses and mixings that span the solution space as shown. 
For sin^2^ greater than ^ 3 x 10~^, z/^ becomes lighter than z/g, and MSW oscillations occur 
between anti- neutrino rather than neutrino states, and are thus not relevant. We note that 

(b) passes through the (small-angle, non-adiabatic) MSW region that is compatible with all 
detectors and also represents a value of m^ (7 eV) that provides a very good fit to COBE 
data in the context of a hot plus cold dark matter scenario. Finally, curve (c) represents a 
singlet mass of 1 TeV and a z/^ mass of ~ 9.8 eV, and terminates where it does because for 
larger mixing angles the z/g becomes heavier than the z/^. Note that the determination of 
which fiavor the z/g oscillates to is made by examining the mixing (diagonalizing) matrix of 
the full ( i.e. tree -|- loop ) mass matrix. A (reasonable) assumption built into the results is 
that z/g is the lightest state. 

We stress that these curves represent a phenomenological exercise more than anything 
else to demonstrate that our calculations can make connection with experiment when the full 
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parameter space, which is quite large, is constrained by imposing physically and empirically 
well-motivated restrictions. 

We note that the singlet mass values chosen by us (100 GeV, 400 GeV and 1 TeV) are 
not in conflict with accelerator |T^ or cosmological |T^ bounds on these particles. 



Finally, we remark that a disparity between the mass scales of the z/g, z/^(z/t-) and that 
of the z/t-(z/^) seems to be required if we take both the solar and COBE implications for 
neutrino masses seriously. In the simple model under consideration here, such a disparity 
arises naturally since the neutrino which contributes to dark matter has a tree level mass 
while the other two have loop masses. 

4 Conclusions 

We have explicitly obtained an analytic form for the radiative two-loop masses acquired by 
doublet neutrinos in models where right-handed singlets are present. We have made an effort 
to keep our calculation general and the expression for the mass correction that we obtain 
may have applications in other models with right-handed neutrinos. We have calculated 
these masses (for the one singlet case) in the light of experimental data from solar neutrino 
detectors and from COBE, within the confines of the MSW solution to the solar deficit. 
By doing so we have made an effort to demonstrate that intermediate scale physics ( i.e. 
physics at < 1 TeV) can lead, in a simple way, to naturally small masses for neutrinos which 
have physically meaningful values, without requiring drastic changes in the presently known 
particle spectrum or gauge group structure. 
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A Appendix: Evaluation of A 



123456 



In this section we discuss the exact evaluation of the fundamental finite two loop four di- 
mensional integral underlying the mechanism. As a first step, though, we consider the more 
general two loop Euclidean space integral, A123456, defined by 

which we will evaluate analytically and then specialize to the case we are concerned with. 
For reasons which we explain below we choose to calculate eqn (A.l) in d-dimensions where 

A-d . 

^ ' (fk (A.2) 



Ik {2tvY 

and /i is an arbitrary mass parameter introduced to ensure the coupling constant remains 
dimensionless in our d- dimensional manipulations. The subscripts on A123456 correspond to 
the masses mf of the integral and we note that the function has certain obvious symmetries, 
Ai23456 = A213456 = A563412, which ought to be preserved in the final expression. The strategy 
to evaluate eqn (A.l) is to use partial fractions to obtain a sum of 2-loop integrals with three 
propagators and then to substitute for the value of each of these sub-integrals, which have 
been considered by other authors in different contexts before, [|1^, [1^, [1^, |T9[. For instance, 
if we define 

'^*^' " XX {p^ + mi){q^ + m']){{p + qf + ml) ^^'^^ 

then eqn (A.l) is built out of a sum of eight such integrals where its only symmetry is Jijk 

= Jjik- Rewriting the numerator of eqn (A. 3) one finds 

Jijk = \[lilj- Ijh- hh- {ml-mj -m])Iijk] (A.4) 

where 

11 



"^■' U, (p2 + m]){{p + qf + m2)(g2 + rnl) ^^'^^ 

and the latter function is totally symmetric, corresponding to a two loop vacuum graph 
(ie zero external momentum). The integral lijk has been considered in |1^, 1^ and a single 
integral representation of it exists, |18, 1^, 2^. For our purposes, however, we have chosen to 



use the elegant formula given in |^ since it is explictly symmetric in the masses. Although 
Ai23456 is itself ultraviolet finite the sub-integrals, eqns (A. 3) and (A. 4), are divergent and 
therefore require regularization. In |19|, |2^ dimensional regularization was introduced to 
control these infinities, which is why we choose to calculate eqn (A.l) in (i- dimensions, so 
that lijk involves double and simple poles in e where d = 4 — 2e. Therefore in the final result 
these must cancel for all mf. As a first step, it is trivial to observe that in the partial fraction 
decomposition of eqn (A.l) the Jj/j type terms, which are also divergent, formally cancel to 
leave only the lijk terms. To proceed we recall the important properties of lijk which have 
been discussed in more detail in |20|. In d-dimensions the exact value, for arbitrary (mass)^, 
X, y and z, is 



i{x,y,z) = i{2a,o,o) + r[Fac-y) + rac 



F{x - \c) 



(A.7) 



where 



^ijk 



T' 



i^r(2 - irf)r(i - irf) 

i[x2 + y^ + z^ - 2xy - 2yz - 2zx\^l^ 
X + y + z 



(A.t 



and 



F{w) 



ds ■ 



a2)(4-'i)/2 



(A.9) 
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The result (A. 7) is valid in the region of {x, y, z) space where a^ > 0. For the case when a^ 
< 0, then the solution is, with b"^ = — a^, 



I{x,y,z) = -/(26,0,0)sin(i7rd) 

+ r[G{lc-x)+G{^c-y) + G{^c-z)] 



where 



and, for example. 



G{w) 



ds ■ 



s2_^52)(4-d)/2 



/(x,0,0) 



T{2 - U)T{3 - d)T\M - l)x- 



d-3 



(A.IO) 



(A.ii: 



(A.12) 



(47r)^r(id)(/i2)'i-4 

which is clearly singular in four dimensions. To obtain the finite part of A123456 each part of 
I{x, y, z) needs to be expanded in powers of e to the 0(1) term and the poles in e cancelled. 
The non-trivial part of this exercise is the e-expansion of the F{w) and G{w) integrals. These 
have been given in [^ and we record that to the e-finite term. 



{ATrfl{x,y,z) 



c _ 1 
2?~7 



3c 



Li 



i[L2-6Li+^(x,y,z) 



+ c(7 + C(2)) + {y + z- x)\ny\nz 

+ {z + X — y)\n.z\n.x + {y + x — z)\n.y\ax\ 



(A.13) 



where C,{n) is the Riemann zeta function, Li = xln x + yln y + zln z, Inx = ln(a;//i^), /i^ 
4:TTe~'^fi'^ and 7 is Euler's constant, and for a^ > 0, 



ax,y,z) = 8a[M(0,) + M(0,) - M(-0,)] 



where 



M{t) = - d(j)\n sinh ( 
Jo 



and the angles (j)x are defined by 



coth"^ 



tC — X 



(A. 14) 



(A.15) 



(A.16) 
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tC — X 



For a^ < 0, then 

where the 9^ angles are given by 

6x = tan~^ 

and L{t) is the Lobachevskij function, 

L{t) = - f rf^lncos^ 
Jo 

Equation (A. 17) can also be rewritten as 

ax,y,z) = 8b[L{e,) + LiOy) - L{-e^)] 



(A. 17) 



(A.18) 



(A.19) 



(A.20) 



where L(t) = J^ dO In cos ^ in order to make the obvious analytic continuation across a? = 
more apparent. It is worth noting that essentially eqn (A.l) has been reduced to a single 
simple function, eqn (A.19), whose properties are well known. We have used the following 
identities in order to write an efficient programme to calculate A123456 for a range of physical 
mass values. For instance, [ETl, 



m 

L{t) 

m 



ivr < t < ivr 



-L{-t) for 

L(i7r - t) + (t - ivr) ln2 - ^L{^7T - 2t) for < t < ivr 

±L(7r±t)^7rln2 



(A.21) 



Therefore, when the argument of the Lobachevskij function is known, the identities of eqn 
(A. 21) mean that one need only write a routine to evaluate L{t) numerically in the range 
[0, ^tt). For example, if < A < 27r then for any integer n 



L(27m + A) = 27rnln2 + L(A) 



(A.22) 



and so on. 
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Returning to the partial fraction form of A123456 with the result for lijk, the c and Li terms 
of the e expansion cancel in the final expression and we can therefore take the limit back to 
four dimensions, e -^ 0. Consequently, we end up with the following analytic expression: 



A 



123456 



4(47r)'^(mf — mD^ml — ml){ml — ml 



irrin 



m^ 



rric 



635 - m^ In ^ In 



mt 



rriq 



TTli. 



-mlln{^,]ln(^ 



mt 



m^ 



- m^ In I ^ I In 






2 1 / "^1 \ 1 
m-^^ m — 2 ^^ 



mt 
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-'"''"(3'"(S-'"^'"(i)'"(Sll <"^^' 

where ^jj^ = C,{mf ,171^1711) and it is evaluated according to eqns (A. 14) or (A. 17) depending 
on whether the particular a^ is positive or negative. A further check on our manipulations to 
obtain eqn (A. 23) is the absence of the arbitrary mass fi which was required at intermediate 
steps to have logarithms whose arguments were dimensionless quantities. 

Although it may appear that the final result is singular in certain cases through denom- 
inator factors like {ml — m\) when m\ = m\, the expression within the square brackets also 
vanishes. Moreover, if one sets m\ = ml + 6, where 6 is small, and expands in powers 
of 6 then in the limit as 5 —^ a non-zero non-singular function of the independent mass 
remains. Further, there is no difficulty with singularities when one or more masses is zero. 
To illustrate this point explicitly we consider the integral A123056 where the zero subscript 
means the corresponding mass of eqn (A.l) is zero. Its form can readily be deduced from 
eqn (A. 23) by taking the m,l -^ limit. However, to do this the behaviour of C,{x, y, z) in the 
z — * limit is required since eqn (A. 23) has terms like \nm\ which are potentially infinite 
in the limit we require. It is is easy to deduce from the explicit representation, eqn (A. 14), 
that 

e(a;, y, z) ^ {x- y) 2U^ (l - f ) + In (^ In (^) (A.24) 

as 2; -^ 0. Thus a little algebra leads to the compact expression, 

1 



A 



123056 



X 



4(47r)''(mf — ml)ml{ml — m^) 

/ 2 2 2\^ / 2 2 2\/- 

/ 2 2 2\/- I / 2 2 2\y- 

(mg - m2 - ^5)^235 + ("^3 - "^2 - "^6)636 

2 _2 _2\ I „l^2 _2 _2\ 



- p(m3,mi,m5)+p(m3,mi,me) 

I / 2 2 2\ / 2 2 2\ 

+ p{m^,m2,m^)- p{m^,m2,mQ) 
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+ A(mi, mg) - X{m^, ttiq) - X{m^, m^) + A(m2, mg) 



(A.25) 



with 



p{x,y,z) = (x-y-z) 



^^"U]^"©^^^"w^n! 



+ zln( - ) In I - 
.xj \y^ 



and 



X{x,y) = (x + y) 



2ix-y)U2{l 



y 



yln 



X 



(A.26) 



(A.27) 



x/ \y^ 

where Li2(t) is the dilogarithm function. Its properties have been discussed extensively in 



22 1 but we make use of the following ones here 



Li2(-t)+Li2(-l/t) 

Li2(t)+Li2(l-t) 



and its integral representation is, [22 



Li2(t) 



-C(2)-iln't fort>0 

C(2)-lntln(l-t) 

rtds 



(A.28) 



/ -ln(l 
JO s 



(A.29) 



where Li2(l) = C(2) = vrV6. 

Finally, another check on our overall expression eqn (A. 23) is the comparison with the 
earlier result of [^ where only ml and ml are non-zero, ie A003400; which was evaluated by 
an independent method. We can easily deduce an expression for A003400 from eqn (A.25) by 
using the relation (A. 24) or by returning to the lijk representation of eqn (A.l) and taking 
the appropriate limits in that case. Useful for the former approach are the properties of the 
dilogarithm function, [^. Whilst in the latter instance we made use of the Tyalor expansion 
of the lijk about zero mass and in particular. 



d^I{x,y,z) 



dydz 



T^{\d-2)T{A-\d)T{b-d)x' 
(47r)^(/i2)'i-4r(irf) 



d-5 



(A.30) 
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whose e expansion is easy to determine. Consequently, we find 

1 frn^X 

^'''''' = "(47r)4(mi- ml) ^"^ 1,4; ^^'^^^ 

This is in total agreement with the explicit calculation of and is a necessary non-trivial 
check that we have the overall normalization of our integral correct, in terms of signs and 
factors of 27r. 
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Figure Captions 



Figure 1: The two-loop diagram which gives rise to the mass corrections considered in this 
paper. 



Figure 2: The MSW solution space for the solar neutrino deficit, from Ref. [|T3|. Superposed 



on it are the 3 curves (a), (b) and (c) which represent sample calculations using our results. 
Each curve shows the mass squared differences and mixings for the two light neutrinos which 
acquire masses radiatively, for fixed values of the masses of the other two neutrinos which 
are massive at tree level. Curve (a) in Figure 2 corresponds to a singlet mass of 100 GeV 
and a Vr mass of ~ 8.6 eV. Curve (b) corresponds to a singlet mass of 400 GeV and a v^ 
mass of 7 eV. Finally, curve (c) represents a singlet mass of ITeV and a v^ mass of ~ 9.8 
eV. 
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This figure "figl-l.png" is available in "png" format from: 



http://arXiv.org/ps/hep-ph/9401329vl 



This figure "figl-2.png" is available in "png" format from: 



http://arXiv.org/ps/hep-ph/9401329vl 



